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It has been shown that the string-flip potential model reproduces most of the bulk 
properties of nuclear matter, with the exception of nuclear binding. Furthermore, 
it was postulated that this model, with the inclusion of the colour-hyperflne in- 
teraction, should produce binding. In some recent work a modified version of 
the string-flip potential model was developed, called the flux-bubble model, which 
would allow for the addition of perturbative QCD interactions. In attempts to 
construct a simple qq nucleon system using the flux-bubble model (which only 
included colour-Coulomb interactions) difficulties arose with trying to construct 
a many-body variational wave function that would take into account the locality 
of the flux-bubble interactions. In this talk we consider a toy system, a mesonic 
molecule in order to understand these difficulties. En route, a new variational 
wave function is proposed that may have a significant enough impact on the old 
string-flip potential model results that the inclusion of perturbative effects may 
not be needed. 



1 Introduction 

For the past 30 years several attempts have been made, with little success, to 
describe nuclear matter in terms of its constituent quarks. The main difficulty 
is the non-perturbative nature of QCD. The only rigorous method for handling 
multi-quark systems to date is lattice QCD. However, this is very computa- 
tionally intensive and given the magnitude of the problem it appears unlikely 
to be useful in the near future. As a result, more phenomenological means 
must be considered. 

A good phenomenological model should be able to reproduce, at least 
qualitatively, all the overall bulk properties of nuclear matter, in particular: 

• nucleon gas at low densities with no power-law van der Waals forces 

• nucleon binding at higher densities 

• nucleon swelling and saturation of nuclear forces with increasing density 

• quark gas at extremely high densities. 

There are many models that attempt to reproduce these properties but none, 
them does so completely. In this talk, only the string-flip potential modeJll'l 
will be considered. This model appears to be promising because it reproduces 
most of the aforementioned properties with the exception of nucleon binding. 

"Speaker. 
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In an attempt to find the root of this problem, a toy mesonic- molecular system 
will be used toiitudy the inclusion of perturbative interactions and variational 
wave functionsQ 

This proceeding is essentially part II of our MRST '94 talkQ For complete- 
ness sake, a brief overview of the previous proceedings will be given, followed 
by our most recent work. 

2 Prom String-Flips to Flux-Bubbles 

In this a section, a brief outline of the generic string-flip potential model is 
presented, followed by an overview of the flux-bubble model. 

2.1 The String-Flip Potential Model 

The string-flip potential model postulates how flux-tubes should form between 
quarks at zero temperature, based on input from lattice QCD and experiment. 
An adiabatic assumption is made, in which the quarks move slowly enough for 
their fields to reconfigure themselves, such that the overall potential energy is 
minimized: i.e., a 



In Eq. |l|, the sum is over all gauge invariant sets, {qm ■ ■ ■ ^n}, of Nq quarks, 
{qi . . . qNg}, in a box of side L with periodic boundary conditions (to simulate 
continuous nuclear matter). 

The physics, or model input, is contained within the many-body terms, 
v(r*rri . . . ?„), in V, which, in general, can be quite complex and virtually im- 
possible to compute. Ironically, the form of v is quite simple: the total minimal 
length, of flux-tubing (or string) connecting a given set of quarks times the 
string tension, a, w ^ ai . Over the past couple of years, many simpliflcations 
of V have been assumed in the form of two-body, three-body, and long com- 
plex chaias-jof quarks, using SUc{2) and SUc{3) with and without "moving- 
colour." floBu All of these models yield similar results for nuclear matter — 

^Typically, a quark is assigned a "fixed-colour" wliicli doesn't change tlirougliout the 
evolution of the gas. Conversely, "moving-colour," is when the quarks are allowed to change 
colour. 




(1) 
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no nuclear binding. A^eeneral synopsis of these models can be found in the 
MRST '94 proceedingsB For simplicity, we shall concentrate on the two-body 
quark-antiquark form, v ^ <^i'qq , with S'[/c(2). 



2.2 The Flux- Bubble Model 



To date, string-flip potential models assumes only non-perturbative interac- 
tions: cf., -v a r ~ ois/r . Flub-bubble models are a recent attempt to include 
perturbative interactions. In particular, they are a derivative of the following 
two-body potentialQ 



ro)0{rij - ro) + OsXij ( — - — ) 



To/ 



0{ro - Tij) , 



(2) 



where Ay is a model dependent colour factor, which is free from Van der Waals 
forces. 

A "flux-bubble" is formed when two or more ends of neighboring flux-tubes 
come within tq of each other, by inserting virtual gg-pairs on the boundary of 
this region, thus leaving a perturbative bubble with non-perturbative flux- 
tubes emanating from it. 

E.g: Consider configuration (a) of quarks, with r > rg, 
about to move to (b), s.t., two of them are within r < rg. 
Then the procedure is to draw a bubble of ro away from 
the two, (b), and to cut the fiux-tubes at the boundary and 
insert virtual qq pairs, (c). Once the potential is computed 
the configuration is restored to (b) before the next move 
is made. 



3) a q 




Virtual qq pair insertion 



Note the insertion of the virtual qq pairs allows the construction of colourless 
objects. These are solely used as a tool to calculate the overall length of the 
flux-tube correctly, and not used in computing the Coulomb term, however, as 
the field energy is already taken into account by the "real" quarks inside the 
bubbles. In general, once the bubbles have been determined, the flvix-tubes 
must be reconfigured in order to minimize the linear part of the potential. 
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Therefore, the most general gq-fiux-bubble potential, with 5J7(2)-colour, 
is given by (c/., Eq. 



ll{qq} l<i 



where Ay equals —3/4 for unlike-quarks and 1/4 for like-quarks. The first 
term represents the string-flip, or linear, potential model, while second term 
contains Coulomb corrections to the this potential (i.e., linear-plus-Coulomb 
potential model) plus new terms describing interactions between the ends of 
the flux-tubes.[] 

2.3 Implementation Problems 

The actual implementation of the flux-bubble potential turned out to be quite 
straightforward. However, trying to find a good many-body variational wave 
function proved to be insurmountable. The main problem was that each trial 
wave function required vast amounts of calculations by hand, writing computer 
code, and lots of computer time (running several machines in tandem) . Human 
CPU time (Thuman ^ 0{l)wk) aside, using a simple mesh minimization scheme 
requires 0{MP/m) of CPU time (tcpu), where M is the mesh size, p is the 
number of parameters, and m is the number of machines. 

A string-flip potential model, for 7 triplets of quarks using fixed SU{2)- 
colour, with a 3-parameter wave function, 

- E (^^^^-)" 

*a/3 = e --{«^"> J]|$(rpJ|, (4) 

Correlation ^ (3 Slater p 

with parameters p (density), f3 (inverse correlation length), and a, would con- 
sume about 5 weeks of CPU time on our 8-nodc farm, here at Carleton, using 
a crude 10"^ mesh. However, because of certain assumptions {a{p) « a{p — 0)) 
and scaling tricks ((/3,p^/'^) — > C(0) (cos 6 , sin 6)) , the CPU time can be re- 
duced to 8hrs. Unfortunately, because of the local nature of the flux-bubble 
model, we cannot use scaling tricks. So, in general, we would require Tcpu 
{10P-^/21)wksl 

Thus, a way of checking different wave functions and minimization schemes 
which do not consume large amounts of CPU time is desirable. In particular, 

'^It is conceivable, that by extending the model to include interactions along the length 
of the flux-tubes, could lead to a model for pion-liquid crystalsH 
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a mini-laboratory is needed in which various aspects of the string-flip and flux- 
bubble potential models, from wave functions to minimization schemes, can be 
investigated, without being too concerned about CPU overhead.'| 

3 Mesonic Molecules 

An interesting place that might make a good mini-laboratory is a mesonic- 
moleculeP Q2, consisting of two heavy-quarks, Q, and two relatively light an- 
tiquarks, q. The heavy-quarks are assumed to be massive enough so as not to 
be disturbed by the motion of the light-antiquarks: i.e., the adiabatic approx- 
imation. By varying the distance, R, between the heavy-quarks a mesonic- 
molecular potential, U{R), can be computed. 

ThCpSchrodinger equation that describes the effective potential, U{R), is 
given bytl 





* = U{R)^ "^^T^ (5) 



where niq = 330 Me^. The potential, V, describes the many-body nature of 
the four quark system and is therefore model dependent. This equation can 
be solved variationally for U{R), at fixed values of R, by guessing the form 
the wave function, and then minimizing U = T + V with respect to the 
parameters in '5. 

Figs. 2. a, shows the results for linear, linear-plus-coulomb, and flux-bubble 
potential models with fixed-colour, using the old string-flip variational wave 
functional 

^^0 = e ^ (7) 

with parameters and (3 . These models yields an average well depth of 
0{3)MeV , which is not enough to bind the Q2 molecule. This suggests that 

''in fact, it was this mini-laboratory wihich enabled us to developed a distributed mini- 
mization algorithm with rcpu ^ 0{p/m)a 

"^Figs. 2. a, also shows results for the semi-relativistic Dirac-Hartree-Fock equation, 

J2 [V? + {m,+<7rp,-)2]vI> = C/(i?)2^ (6) 

min{Qq} 

(c/., |), for m, = and 330 MeV. 

•^Analysis was also done for fixed values of a , at 2 and 1.74, which shows a 1% variation 
in well depth, suggesting past assumptions about a are valid. 
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none of these models will lead to nuclear binding^ or that perhaps something 
is wrong with the wave function itself. 




Fig. 2: Binding energy as function of heavy-quark separation for the old (a) and new 
(b) wave functions with fixed-colour, and respective plots (c) and (d) for the flux-bubble 
model with moving-colour. Shown are the string-flip (linear and linear-plus-Coulomb) 
and flux-bubble (plus an SUc{3) composite model) models for both wave functions, and 
semi-relativistic plots for the old wave function with flxed-colour. 



Although the Q2 system is far removed from its H2, hydrogen molecule, 
cousin from a dynamical point of view and the motivations for achieving lo£ 
ization are quite different, it would seem plausible to use a similar ansatzt 

^a,p^e +e , (8) 



^We are not surprised that the flux-bubble model yields very littlipichange, as the it was 
postulated that the colour-hyperfine interactions are more importantu 
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where setting a = 1 yields the H2 case. Now, if qiQi and 52 (92 represent 
two separate mesons, say, then the first term represents the internal meson 
interactions while the second term represents the external meson interactions. 
Notice that the external interactions shut off as the separation, R , between 
the two heavy-quarks becomes large, which is the desired property. 

Fig. l.b, shows the results for linear, linear-plus-coulomb, and flux-bubble 
potential models with fixed-colour, using the new _ff2-like variational wave 
function]] The contrast between the results for the old and new wave functions 
is quite dramatic! The new wave function gives a factor of 0(27) increase in 
well depth; deep enough to bind a heavy-quark system with a reduced-mass 
greater than O{660)Mev)i 

Figs, l.c and l.d show the results for the flub-bubble models with moving- 
colour for the old and new wave function respectively. Fig. l.c shows a lin- 
ear core (with slope a) with a repulsive barrier followed by an abrupt cut-off 
at 1.5/m. This suggests that minimal energy state is a system consisting 
of two mesons, one containing heavy-quarks and the other containing light- 
antiquarks. Fig. l.d shows potential with a very strong Coulomb-like (i.e., 
a ~ 0(125.5)71^6^/771) interior, about SOOMeV deep, with a Yukawa-like ex- 
terior. This suggests a very tightly bound Q2 system, however, the Q2-well 
for the old wave function is much deeper, suggesting the aforementioned two- 
meson system is the minimal energy state. Of course this is not physical, and 
suggests that phenomenological models using SUcC^), with moving-colour, will 
yield unsatisfactory results. 

Fig. l.b, shows the results for an SUd^) composite model, a "baryonic 
molecule," which assumes the two heavy-quarks are actually a composite of 
two light-quarks (also, the orbiting antiquarks have been replaced with quarks). 
This model avoids dissociation into mesons, since a flux-tube can no longer be 
formed between the (composite) heavy-quarks, thus, restoring the potential to 
its generic (original) form. 

It should be stressed that the Q2 system is a toy-model, a pedagogical aid, 
used to gain better insight into the complex nature of nuclear interactions. 
Nuclear systems (at low temperature) consist of equal mass quarks, whereas, 
the Q2 consists of unequal mass quarks. Therefore, what is applicable for Q2 
system may not translate very well to nuclear matter. That said, the results 
herein, suggest that the many-body wave function for nuclear matter should 
be of the form. 



''Also shown, are results for a fixed at 2 and 1.74 . 



^' Perm 
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which is a product of a totaUy symmetric correlation wave function times a 
totaUy antisymmetric Slater wave function. A simple test of this, ansatz, would 
be to consider in the context of a qq gas with fixed- S'[/c (2). 

4 Conclusions 

The Q2 system has proven to be a very useful aid for trying to sort out the 
complexities of model building for nuclear matter. The details of the mechan- 
ics, from wave functions to dynamics to practical computing methods, of the 
flux-bubble model have been throughly investigated. However, it remains to 
be see how applicable this work will be to nuclear matter; something that we 
hope to investigate in the near future. 

In closing, it appears that the flux-bubble model may prove not only useful 
for modeling nuclear matter, but also useful for modeling mesonic molecules 
as well. 
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